Spin-Peierls transition of the first order in S = 1 antiferromagnetic Heisenberg chains 
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We investigate a one-dimensional S — 1 antiferromagnetic Heisenberg model coupled to a lattice 
distortion by a quantum Monte Carlo method. Investigating the ground-state energy of the static 
bond-alternating chain, we find that the instability to a dimerized chain depends on the value of the 
spin-phonon coupling, unlike the case of S = i. The spin state is the dimer state or the uniform 
Haldane state depending on whether the lattice distorts or not, respectively. At an intermediate 
value of the spin-phonon coupling, we find the first-order transition between the two states. We also 
find the coexistence of the two states. 
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I. INTRODUCTION 

A spin-Peierls system has been one of the most fasci- 
nating topics of the low-dimensional quantum spin sys- 
tem. The one-dimensional S = g antiferromagnetic 
Heisenberg chain shows a spontaneous lattice dimeriza- 
tion due to the spin-phonon coupling at low tempera- 
ture, where the strength of the antiferromagnetic cou- 
pling alternates between strong and weak values. For 
S — 5, there have been numerous studies on the ground 
state of the spin-Peierls system. It has been shown by 
a bosonization method that the bond alternation 
Ji = J[l + ( — l) l S] causes the magnetic energy gain pro- 
portional to <5 4 / 3 for small S. Therefore, the magnetic 
energy gain always exceeds the elastic energy for small 
6, because the elastic energy is proportional to S 2 . Thus, 
the dimerized state is always realized in the ground state 
regardless of the value of the spin-phonon coupling. This 
mechanism has been confirmed in numerical works ||,|J , 
where the dimerized lattice configuration is obtained in 
the balance of the magnetic energy and the elastic en- 
ergy in the ground state. Properties at finite temperature 
have been also investigated |j-[7]] , where the distribution 
of the lattice distortion in the thermal equilibrium state 
is obtained. Recently, impurity-induced long range or- 
ders have been found and their physical origin has been 
discussed extensively ||. 

On the other hand, only a small amount of study on 
the spin-Peierls transition for S > h cases has been done. 
In the case of integer spin antiferromagnetic Heisenberg 
uniform chains, as pointed out by Haldane (9J, an energy 
gap exists between the ground state and the first ex- 
cited state, and the spin-spin correlation function decays 
exponentially contrary to the case of half-odd-integer 
spin chains. Because of the qualitative difference of the 
ground-state property between half-odd-integer spin and 
integer spin, the nature of the spin-Peierls transition in 



integer spin systems is expected to be different from 
that in half-odd-integer spin systems. Guo et al. in- 
vestigated the magnetic energy gain due to the dimer- 
izaion for S = 5, 1, § and 2 by an exact diagonahzation 
method JIo[| . They pointed out that the transition in sys- 
tems with integer spin occurs only for large value of the 
spin-phonon coupling while the transition occurs regard- 
less of the value of the spin-phonon coupling in systems 
with half-odd-integer spin. There this difference was dis- 
cussed from a view point of the valence-bond-solid (VBS) 
picture. However, the nature of the phase transition was 
not discussed in detail. 

In the spin-5 antiferromagnetic Heisenberg bond- 
alternating chains, successive quantum phase transitions 
occur at 25 gapless points in the range — 1 < 5 < 
1 [O-M]. The transitions can be considered as ones be- 
tween the different VBS states, as shown in Fig. [I]. In 
the case of S — 1, the continuous transition between the 
Haldane state and the dimer state occurs at the critical 
point S c ~ 0.26 [||. 

In this paper we study the case of S = 1 by a quantum 
Monte Carlo method. First we consider the instability to 
a dimerized chain in the ground state. In the static bond- 
alternating chain, the ground state depends on the am- 
plitude of the bond alternation as we mentioned above. 
Therefore, there are two possibilities for the instability 
to the dimerized chain, namely, the chain is unstable to- 
ward the dimerized chain with a small bond alternation 
where the spin state is the Haldane state, or toward the 
dimerized chain with a large bond alternation where the 
spin state is the dimer state. To which state the system 
changes depends on the value of the spin-phonon cou- 
pling. We also find that the uniform Haldane state and 
the dimerized state coexist in the ground state at an in- 
termediate value of the spin-phonon coupling. On the 
other hand, we find that the dimerized state is more sta- 
ble than the uniform Haldane state at finite temperature, 
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which can be attributed to an entropy effect. Further- 
more, we investigate an open chain with fixed boundary 
bonds, by which we manipulate the configuration of the 
bonds at edges. Locating bond alternation at one end 
and uniform configuration at the other end, we force the 
system to have the two regions, and observe the situa- 
tion where the uniform Haldane state and the dimerized 
state coexist in the chain. We investigate the domain 
wall between the two regions. 

This paper is organized as follows. In Sec. [n] we explain 
a model and a numerical method. We also refer to the or- 



der parameters. In Sec. Ill the magnetic energy gain due 
to the dimerization is investigated in order to study the 
instability to a dimerized bond structure in the ground 
state. In Sec. the distribution of the bond distortion 
at finite temperature and the effect of the entropy on the 
coexistence of the two states are investigated. In Sec. |v| 
the domain- wall configurations between the uniform Hal- 
dane state and the dimerized state are investigated. We 



summarize our results in Sec. VI 
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FIG. 1. The conceptual configurations of the VBS state 
for 5=1. 



II. MODEL AND METHOD 

The spin-Peierls system is a spin-phonon coupled sys- 
tem described by 
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where the distortion of the exchange coupling is as- 
sumed to be proportional to the lattice distortion, that is, 
AJi/J — aAi. Here a is a coefficient that gives the dis- 
tortion of the exchange coupling and is negative (a < 0) 
because the shrink of the lattice causes the increase in 
the exchange coupling. In this paper, for the simplicity 
of notation, we use the change of the exchange coupling 
Aj = aAi instead of A 4 as a variable for the distortion 
and also scale the elastic constant as k = k/a 2 , which 
leads to 
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We study this Hamiltonian (g) for S — 1 by a world 
line quantum Monte Carlo (QMC) method in order to 



investigate thermodynamic properties of the model. We 
use a quantum Monte Carlo loop algorithm with contin- 
uous imaginary time [l8| . This algorithm is effective 
for studying the property at very low temperature. We 
have studied the S = I system with the bond distor- 
tion by the loop algorithm, where we investigated how 
the dimerization develops as a function of the tempera- 
ture and also the effect of the bond fluctuation at high 
temperature where the averaged bond structure is still 
uniform jjj . The system is updated in the following way. 
Starting from an arbitrary spin and bond configuration, 
the spin configuration is updated by the loop algorithm 
with the fixed bond configuration. The loop algorithm 
of a general spin is constructed as the extension of the 
S = h format The bond configuration {A,} is up- 
dated by the Metropolis algorithm with the fixed spin 
configuration. Here we treat the bond distortion as a 
classical quantity, although the quantum phonon effect 
is also an interesting problem f|[|. These two updated 
procedures are carried out alternately, therefore we ob- 
tain the thermal equilibrium state both of spin and bond. 
In the present study, we impose the following constraint 
in the simulation, 
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in order to fix the total length of the lattice. Besides, we 
restrict the bond distortion to — 1 < A, < 1 in order to 
confine the exchange coupling to be antiferromagnetic. 
The periodic boundary condition is adopted except in 
Sec. where we adopt the open boundary condition 
with fixed boundary bonds. We set the Boltzmann factor 
fee = 1 and use it as the unit of the energy. We fix the 
uniform exchange coupling J = 1 and investigate cases 
with various values of k. 

We investigate the ground state energy of static bond- 
alternating chains in order to determine how much en- 
ergy the system gains from the magnetic interaction in a 
given bond configuration. In order to obtain the energy, 
we perform a QMC method without a bond update. The 
initial 10 3 Monte Carlo steps (MCS) are discarded to ob- 
tain the thermal equilibrium state and the data of physi- 
cal quantities are sampled in the following 10 4 MCS. For 
the study of bond- fluctuating chains, larger simulations 
such as 10 6 MCS are performed to obtain good conver- 
gence of the data. The sampled data is divided into 10 
bins, and the error bar is estimated from the standard 
deviation in this set of data. 



A. The order parameters 

In this paper we study orderings of the bond configu- 
ration and the spin configuration. In order to investigate 
the bond ordering, we introduce the bond staggered order 
parameter 
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which represents how the bond alternates in the whole 
chain. When the bond configuration takes an alternat- 
ing arrangement on average, the spin state would change 
according to the amplitude of the bond alternation. In 
the static bond-alternating chain, the hidden Z2 x Z2 or- 
der exists in the Haldane state, while not in the dimer 
state. The hidden order is detected by the string order 
parameter introduced by den Nijs and Rommelse [f20[, 



Ostring = lim C st ring(«, j), 
|»— j|— ►oo 



where 



C st ring(i, j) = S* exp 
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In order to investigate whether the spin state takes the 
Haldane state or the dimer state, we observe the string 
correlation by the following quantity (2TJ , 
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where N s is the number of nonzero spins, and N s 
£*Li \Sf\- Then we calculate (0£ R ): 
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which gives an estimate of the string correlation length 

^string- 



monotonously as 5 increases. Besides, we find a singular 
behavior in E sp i n (S) at the critical point 5 C = 0.26 where 
the phase transition occurs between the Haldane state 
and the dimer state |T^| . The magnetic energy curve can 
be fitted well by an even function in each region of 5, 



E spin {S)/N = 



bS 2 (5 < S c ) 
-b'5 2 + c'5 4 (5>6 C ), 



(10) 



with (a, b) = (-1.40, -0.385) and (a', 6', c') = 
(-1.38, -0.770, 0.148). Adding the elastic energy, 
^lattice (8)/N = ^k5 2 , we obtain the ground-state energy 
as a function of 5, 



E s (6) = E spin (5) + £lattice(<5), 



(11) 



which is interpreted as an effective potential for the bond 
distortion. The results are shown for some values of A; in 
Fig. ||. First let us consider the case of small S. When 
hk > \b\, the uniform configuration 5 = gives a stable 
point at least locally. On the other hand, when ifc < 
the uniform configuration becomes unstable. Here the in- 
stability to the dimerized chain occurs at |fe = |fe|. How- 
ever, if we look at the shape of E s (6) globally, we find 
that a minimum point of the potential jumps from 5 = 
to 5q (> S c ) discontinuously. Namely, the energy shows 
the characteristic behavior of the first-order transition. 
The fourth term c'<5 4 suppresses the divergence and gives 
a finite value So, which minimizes the total energy. In 
the present study, however, we restrict the value of 8 to 
\5\ < 1. Therefore, 5=1 gives the minimum value of the 
total energy when <5 > 1. This constraint is introduced 
in order to confine the exchange coupling to be antifer- 
romagnetic. Physically, this constraint corresponds to a 
nonlinear increase of the energy of the lattice distortion 
or a nonlinear relation between the lattice distortion and 
the change of the exchange coupling. 



III. DIMERIZED GROUND STATE 

In this section we investigate the dependence of the 
ground-state energy on the amplitude of the bond dis- 
tortion by the QMC method. This dependence gives an 
effective (or adiabatic) potential for the bond distortion, 
by which we can study the instability to the dimerized 
chain in the ground state. Here we consider static dimer- 
ized chains with an amplitude 5, 

A, = (-iyS. (9) 

In Fig. |[ we show the magnetic energy as a function 
of 5, E sp i n (5) for the system of N = 64 and T = 
0.01. By studying the size and temperature dependen- 
cies, we confirm that this system size is sufficiently large 
enough to obtain the ground-state energy without the 
finite-size effect, and this temperaure is sufficiently low 
enough to obtain the ground-state energy without the 
finite-temperature effect. The magnetic energy decreases 
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FIG. 2. The 5 dependence of the ground-state magnetic 
energy with the dimerization Aj = ( — l) l <5. As a guide to the 
eye, the values for 5 < 5 C and for 5 > S c are plotted with 
the solid circles and the open circles, respectively. The solid 
lines denote fitting curves for each region of 5 on both sides 
of 5 C = 0.26 as mentioned in the text. 
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FIG. 3. The 5 dependence of the ground-state energy for 
various values of k. 
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FIG. 4. The k dependence of So, which minimizes the 
ground-state energy. The solid line denotes a stable value as 
a function of k. The dotted line denotes an unstable value as 
a function of k. Here we restrict 8 to within S < 1. 

The ground-state value 5q as a function of k is plot- 
ted in Fig. |4. For large values of k, the energy increases 
monotonously as 5 increases and the ground state is the 
uniform Haldane state. For intermediate values of k, 
E g (S) has two minimum points. The minimum 5 = 
corresponds to a uniform chain where the spin state is 
the Haldane state, while <5 = <5o (> S c ) corresponds to a 
bond-alternating chain where the spin state is the dimcr 
state. In particular the energies of the two minima are 
the same at k = k c — 1.30. Here two spin-lattice con- 
figurations coexist in the ground state. For small values 
of k, the energy decreases monotonously and the ground 
state is the dimerized state. It should be noted that <5o 
is always larger than S c . Therefore, we conclude that 
the spin state is the dimer state whenever the latice is 
spontaneously distorted. Namely, in the ground state of 
the present system, the Haldane state appears only in 
a uniform chain and does not appear in a weakly bond- 
alternating chain. 

In the coexistent range, the maximum of the energy be- 
tween the two minima corresponds to an energy barrier 
AE between the uniform Haldane state and the dimer- 
ized state. In particular around k = 1.3, the energy bar- 



rier is rather small, i.e., AE ~ 0.02. Therefore, a large 
fluctuation is expected at finite temperature higher than 
the height of this energy barrier, i.e., T ~ 0.02. 



IV. THERMAL FLUCTUATION 

In this section we study the bond distortion {Ai} at fi- 
nite temperature by the QMC method. First let us study 
the distribution of the bond distortion. At T = 0, the 
distribution P(A) is the delta function with one peak at 
A = for the uniform chain or two peaks at A = ±<5o 
with the same amplitude for the bond-alternating chain. 
The peaks are broadened at finite temperaure due to 
the thermal fluctuation. At high temperature, we find 
a broad single peak near at A — 0. On the other hand, 
at low temperature, we find a sharp distribution with one 
peak at A = or two peaks at A = ±<5o according to the 
value of k. Thus we confirm that the dimerized ground 
state obtained in the previous section is actually realized, 
and it ensures that the lattice takes the dimerized config- 
uration although many other configurations are possible 
in principle. 

As mentioned in the previous section, the uniform Hal- 
dane state and the dimerized state coexist in the ground 
state at k = k c . It is expected that the distribution is 
given by a combination of each distribution. However, 
the system takes either state at a time in the simulation. 
In Fig. ^ we show the distribution function for k = 1.3 
at a very low temperature T = 0.01. We find that the 
system is trapped either in the uniform state [Fig. ||(a)] 
or in the dimerized state [Fig. ||(b)]. Once the system is 
trapped in one state, it is rather difficult to transfer to the 
other state at the temperature lower than the energy bar- 
rier. In Fig. ^|(b), we find an asymmetrical distribution 
for the dimerized state. The amplitude at the positive 
A peak is larger than that at the negative A peak. This 
difference can be attributed to the energy difference of 
the state. We have observed this asymmetry in the case 
of S = \ as well §. 

On the other hand, when the temperature exceeds the 
energy barrier, a change of states between the uniform 
state and the dimerized state occurs due to the thermal 
fluctuation. In Fig. we show the time evolution of 
the bond staggered order Ag g and the string order 0£ R 
for k = 1.3 at a modestly low temperature T = 0.03. Al- 
though the coexistence of the two states is expected in the 
ground state at k — 1.3, the system stays in the dimerized 
state through almost the whole simulation. Namely, the 
dimerized state is favored at a finite temperature. This 
must be due to an entropy effect. In order to make it 
clear, we study the case of k = 1.33. As shown in Fig. [?|, 
we find jumps between the two stable states. Namely, the 
bond configuration takes a uniform configuration and an 
alternating arrangement by taking turns in the simula- 
tion, and correspondingly the spin state takes the Hal- 
dane state and the dimer state. At this value of fc, the 
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ground-state energy in the dimerized state is larger than 
that in the uniform Haldane state. However, the system 
jumps between the two states almost equally. This in- 
dicates the same free energy of the two states, which is 
brought by the exquisite balance between the quantum 
fluctuation in the ground state and the thermal fluctua- 
tion at finite temperature. We indeed find a distribution 
with three peaks as shown in Fig. |^. 
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FIG. 5. The distribution function of the bond distortion 
at T — 0.01 for k = 1.3 and N = 32. (a) corresponds to the 
uniform Haldane state and (b) corresponds to the dimerized 
state. 
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N=32, k=1.3, T=0.03 : 
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FIG. 6. The MCS dependence of the bond staggered order 
A 2 g and the string order 0£ R at T — 0.03 for k = 1.3 and 
N = 32. The whole simulation is divided in blocks of 10 4 
MCS and the thermal averaged value is calculated in each 
bin. 



This entropy effect can be understood from the shape 
of the effective potential E g (S) shown in Fig. ||. If we 
approximate the potential shape around the bottom of 
the minima, we have effective stiffness constants an and 
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where AE(S ) = [E g (S ) - E K (0)} /N. From Eq. QIC 
we find an = 0.560 and au = 0.480 for the critical value 
k c = 1.30. The free energy of the system in the harmonic 
approximation is given by 
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Thus the difference of the free energy is 
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FIG. 7. The MCS dependence of the bond staggered order 
and the string order Ol R at T = 0.03 for k — 1.33 and 
N — 32. The whole simulation is divided in blocks of 10 4 
MCS and the thermal averaged value is calculated in each 
bin. 
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FIG. 8. The distribution function of the bond distortion 
at T = 0.03 for k = 1.33 and N = 32. 
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If we set the value of k to be equal to the criti- 
cal value k c = 1.30, then AE(8 ) = and AF = 
\kBT\n(aH I 'od) > 0. Therefore, the dimerized state 
is more stable than the uniform Haldane state at finite 
temperature. In order to compensate for the entropy 
term, we can set Ed(8o) > 0. Indeed, if we use k — 1.33, 
the difference of the free energy becomes very small at 
T = 0.03, which is consistent with the above observation. 



V. COEXISTENT STATE 

In this section we study the domain- wall structure be- 
tween the uniform Haldane state and the dimerized state. 
As we mentioned in the previous section, in the periodic 
chain without the inhomogeneity, the bond configuration 
takes the uniform configuration or the dimerized config- 
uration at low temperature, and it is rather difficult to 
see the change between the two states. In order to see 
the domain- wall structure, we fix different bond config- 
urations at each chain end: a uniform configuration in 
the left side and an alternating configuration in the right 
side. The bonds between them fluctuate and take a ther- 
mal equilibrium configuration. In order to investigate 
spin and bond configurations in the coexistent state, the 
amplitude of the bond alternation of the fixed bonds is 
set to the amplitude of the dimerized state that corre- 
sponds to the minimum in the ground-state energy curve 
obtained in Sec. 
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In Fig. |9|, we show the bond configuration at a very 
low temperature T = 0.01 for k = 1.3. We find the 
domain wall between the uniform and the dimerized con- 
figurations. Namely, the bonds relax into the uniform 
configuration in the left side and into the dimerized con- 
figuration in the right side. In the uniform region, the 
lattice shrinks uniformly where the strength of the ex- 
change coupling is larger than that of the original uni- 
form chain. This shift of the exchange coupling causes 
the magnetic energy gain. On the other hand, in the 
dimerized region, A, changes alternately. However, the 
amplitude of negative A increases, which compensates 
for the increase in the uniform region. According to the 
bond configuration, the spin state takes a characteristic 
state in each region. In Fig. [lO[ we show the string cor- 
relation function from the left edge site. In the region of 
the uniform bond configuration, the spin state takes the 
Haldane state and the string correlation persists. On the 
other hand, in the region of the dimerized bond configu- 
ration, the spin state takes the dimer state and the string 
correlation decays exponentially. Thus we find that the 
uniform Haldane state and the dimerized state coexist 
in the chain. The domain wall between the two states 
shows a soliton structure. Finally let us study the local 
magnetic response of the chain. In Fig. 
local-field susceptibility, 
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FIG. 9. The bond configuration in the open chain with 
fixed boundary bonds at T = 0.01 for k — 1.3 and N = 64. 
The fixed bonds are enclosed in circles. 
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FIG. 10. The string correlation function from the left edge 
site in the open chain with fixed boundary bonds at T = 0.01 
for k — 1.3 and TV = 64. Fixed boundary bonds are the same 
as that in Fig. []. 
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FIG. 11. The local-field susceptibility x'° cal in the open 
chain with fixed boundary bonds at T = 0.01 for k — 1.3 
and N = 64. Fixed boundary bonds are the same as that in 
Fig. | 

which represents the degree of quantum fluctuation at 
each site. We find that the response is strong at edges 
of the region of the Haldane state, which is similar to a 
pure open Haldane chain p2[. 
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VI. SUMMARY 

In this paper, we investigated the instability to a 
dimerized chain of a one-dimensional S = 1 antiferro- 
magnetic Heisenberg model coupled to a lattice distor- 
tion by a quantum Monte Carlo method. From the S 
dependence of the ground-state energy E s (8), we found 
that the instability to the dimerized chain depends on 
the value of the spin-phonon coupling unlike the case of 
S = |. We also found the first-order transition between 
the uniform Haldane state and the dimerized state. We 
investigated spin and bond configurations in the coexis- 
tent range. The system takes turns jumping between the 
two states in the simulation. By studying the tempera- 
ture dependence of the value of k at the coexistence, we 
found that the dimerized state is favored at finite tem- 
perature due to an entropy effect that can be explained 
in a harmonic approximation of the effective potential 
function E(5). We also studied the domain- wall struc- 
ture between the two states in an open chain fixing the 
edge bonds to have different states. 

As mentioned above, we found a peculiar character in 
the case of S = 1 that is not observed in the case of 
5=5. In a realistic material, this phenomena could be 
observed as the spin-Peierls transition for small values of 
k, namely, large values of the spin-phonon coupling. We 
also expect that the domain wall between the Haldane 
and the dimer spin states will be found in some material. 
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